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We address the role of fluctuations in strongly interacting matter during the dense stages of a heavy- 
ion collision through its electromagnetic emission. Fluctuations of isospin charge are considered in 
a thermal system at rest as well as in a moving hadronic fluid at fixed proper time within a finite 
bin of pseudo-rapidity. In the former case, we use general thermodynamic relations to establish a 
connection between fluctuations and the space- like screening limit of the retarded photon self-energy, 
which directly relates to the emissivities of dileptons and photons. Effects of hadronic interactions 
are highlighted through two illustrative calculations. In the latter case, we show that a finite 
time scale r inherent in the evolution of a heavy-ion collision implies that equilibrium fluctuations 
involve both space-like and time-like components of the photon self-energy in the system. Our study 
of non-thermal effects, explored here through a stochastic treatment, shows that an early and large 
fluctuation in isospin survives only if it is accompanied by a large temperature fluctuation at freeze- 
out, an unlikely scenario in hadronic phases with large heat capacity. We point out prospects for the 
future which include: (1) A determination of the Debye mass of the system at the dilute freeze-out 
stage of a heavy-ion collision, and (2) A delineation of the role of charge fluctuations during the 
dense stages of the collision through a study of electromagnetic emissivities. 

PACS: 12.38.Mh, 25.75.-q, 13.40.-f, 24.10. Nz 



I. INTRODUCTION 

Charge fluctuations extracted through event-by-event 
measurements in high-energy heavy-ion collisions have 
recently attracted attention as a potential signature of 
Quark-Gluon Plasma (QGP) formation [Q. In the con- 
fined phase charges are carried by hadrons in integer units 
of e, while in the deconfmed phase they are attached 
to quarks in fractional units. This leads to charge fluc- 
tuations and associated thermodynamic susceptibilities 
that are distinctly different in the two phases 0. If sig- 
nificant primordial fluctuations survive the evolution of 
the collision, one might hope to deduce the properties of 
deconfmed matter at the early stages, analogous to the 
manner in which primordial fluctuations in the early uni- 
verse are inferred through the detection of anisotropy in 
the cosmic background radiation. Beginning with the ex- 
pected signals in rapidity distributions || , many sources 
and signals of fluctuations in hadronic observables have 
been considered (cf. [Q and references therein). For more 
recent discussions, see, for example, Refs. pf||- 

In experiments at the CERN SPS and RHIC to date, 
observed fluctuations in hadronic multiplicities, kaon to 
pion ratios, and transverse momenta agree reasonably 
well with model calculations of these quantities at the di- 
lute freeze-out stage of the system 0. Such studies, how- 
ever, have not unraveled the physics associated with the 
QGP phase transition. Whether future collider experi- 
ments at RHIC and LHC will reveal tell-tale signatures 
in fluctuations of hadronic observables alone remains to 
be seen. 



It must be noted that the measured fluctuations in 
hadronic observables result from a combination of sev- 
eral factors that include: (1) The evolution of the created 
system from a dense partonic stage to a dilute hadronic 
stage. The duration of each stage is especially criti- 
cal inasmuch as the memory of the partonic stage is 
more likely to be retained by the hadrons provided the 
hadronic stage is of relatively short duration; and (2) 
The extent to which interactions are different in the two 
phases. A caveat, of course, is that if duality between the 
partonic and hadronic worlds holds, equivalent descrip- 
tions of observables may be achieved by considering inter- 
actions in the two phases. A good example is provided by 
the cross sections ratio R = <r e+e - ^ hadrons /<r e+e - /1 - 
for yfs > 1.5 GeV; in this case, interacting hadronic 
states yield a description that is equivalent to one that 
involves perturbative quarks and gluons. 

Our principal objective in this paper is to empha- 
size that possibly the only way to address the role of 
fluctuations during the dense stages of the collision is 
to exploit electromagnetic probes in combination with 
hadronic probes. Towards this goal, we study relations 
between electromagnetic radiation and thermal fluctua- 
tions of isospin charge in a system at rest as well as in 
an expanding hadronic fluid at fixed proper time and in 
appropriately chosen bins of pseudo-rapidity. 

For a system at rest, we utilize general thermody- 
namic relations to establish connections between fluctu- 
ations in matter and the space-like screening limit of the 
retarded photon self-energy, a quantity which directly 
controls the thermal emissivities of dileptons and pho- 
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tons. We consider two examples that exemplify the ef- 
fects of hadronic interactions on electromagnetic radia- 
tion. In addition, we point out the close connections 
that exist between fluctuations, susceptibilities, and aver- 
age phase space densities. At the dilute freeze-out stage, 
the average phase space densities of pions and kaons can 
be inferred through intcrfcrometric (Hanbury-Brown and 
Twiss) studies |9|,|l0|] and provide for a direct determina- 
tion of the system's Debye mass. 

For a system undergoing relativistic expansion, we de- 
rive a general formula for fluctuations in isospin charge 
binned in pseudo-rapidity intervals of width A. The con- 
stituents of the moving fluid, be they partons or hadrons, 
are assumed to be in local thermal equilibrium. In the 
case of hadrons, we include finite pion and baryon chemi- 
cal potentials. We cast our result in a readily usable form 
for the analysis of experiments, which have unavoidable 
restrictions in rapidity and transverse momenta due to 
limited acceptance in these variables. We show that finite 
time scales t inherent in heavy-ion collisions imply that 
equilibrium fluctuations involve both space-like and time- 
like components of the photon's polarization function 
with predominant contributions arising from frequency 
/ sa 2it/t and wavelength A w 7rrsh(A/2). We also 
consider effects of non-equilibrium (for previous work, 
see for example, Ref. fll ) with a view to explore the 
consequences of possible temperature fluctuations at the 
freeze-out stage. 

This article is organized as follows. In Sec. y, the 
relationship between isospin fluctuations and space-like 
physics that leads to screening of charges is established 
in an infinite volume system at rest and in thermal equi- 
librium. Sees. [n|.A through ^.D briefly summarize known 
relations that have enabled us to interrelate fluctuations, 
susceptibilities, and average phase space densities. In 
Sec. ||.E, connections to electromagnetic (dilepton and 
photon) emissivities are set up. Two illustrative calcu- 
lations incorporating effects of hadronic interactions in 
the dilepton and photo n emission rates are presented in 
Sec. |n].F. In Sec. Ill A. the role of fluctuations in heavy- 



ion collisions, in which the system created expands from 
an initially dense partonic system to a dilute hadronic 
system prior to freeze-out, is discussed with particular 
emphasis on electromagnetic probes. We demonstrate 
that in a relativistically expanding fluid, thermal fluc- 
tuations probe both spac e- like (screening) and time-like 
(emission) physics (Sec. Ill A). A discussion that re- 
lates o ur find ings to those in earlier works is contained 
in Sec. 
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In Sec. [V 



we show how isospin fluctu- 
ations present in the early stage can only survive if a 
sizable temperature fluctuation is simultaneously present 
at freeze-out. This section contains an analysis of non- 
equilibrium effects treated stochastically through addi- 
tive Gaussian and multiplicative power-law heat flow. 
Our summary and conclusions are contained in Sec. [v|. 



II. THERMALLY EQUILIBRATED SYSTEM AT 
REST 

A. Fluctuations 

We begin by considering a system of infinite spatial 
3-volume V3 at rest and in thermal equilibrium. The 
isospin charge operator, Q a , is obtained from the zeroth 
component of the pertinent current, V°, through 



Q Q (0) 



d 3 ccVo(0,f) 



d 3 x q{x)-io — q{x) , (1) 



where q denotes quark-field operators and f are the stan- 
dard SU(2) isospin matrices (here and in what follows 
calligraphic letters indicate operators). The thermal fluc- 
tuations are defined by 

AQ 2 = (Q a (0) Q a (0)) c 

= { d 3 xd 3 x' (VS(0,x)Vg(0,x')) c 



= V 3 / d 3 x<V a (0,a J )Vo°(0,0)) c 



= V 3 ~i I d 3 x I dx°9(x°) ([Vg(x),VS(0)]) , (2) 



where 



(O) = Tr 



e -P(H-iJ,Q-n)Q 



(3) 



indicates the thermal expectation value of the operator 
O. The subscript c in (• • •) refers to connected parts. The 
relevant chemical potential is denoted by fx, j3 = 1/T is 
the inverse temperature, and f2 is the thermodynamic 
potential. The last line in Eq. (0) follows from general 
analytic properties of the retarded correlator with the 
order of integration fixed to enforce approach from the 
space-like regime. 



B. Susceptibilities 

The local fluctuations in a globally conserved quan- 
tity can be expressed through a derivative of the thermal 
expectation value of O with respect to the associated 
chemical potential (see, e.g., Refs. p2|JI^]): 



(AO) 2 = (O 2 ) - {Of 



-T 



d(Q) 
<9/i 



(4) 



Since the charge-density is defined as a partial derivative 
of the thermodynamic potential, its thermal fluctuations 
are given by the charge susceptibility 



Xch 



1 d 2 fl _ (AQ) 2 



v 3 dtf 



ch 



TV 3 



(5) 
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A general relation between the susceptibility and the 
00-component of the electromagnetic (e.m.) polarization 
tensor is (see, e.g., Refs. [|Lp#jl(|) 



Xch 



n°°{q =0,<f^0) 



(6) 



where mo is the Debye mass (strictly speaking, this re- 
lation is valid only to leading order in e 2 ). Equation (^), 
also implied by Eq. (|J), then provides a direct link to 
the electric screening length in the system. On kinemati- 
cal grounds, one expects m D (T) to receive contributions 
chiefly from low-lying mesonic excitations in a hadronic 
gas and from (nearly) massless quarks and anti-quarks in 
a partonic gas. 



C. Examples in Hadronic and Quark-Gluon-Plasma 
Phases 



MeV. As pointed out in Ref. H, inclusion of higher- 
mass hadronic resonances reduces this number further 
to R % h G = 0.06 — 0.07. This is readily understood from 
the fact that heavier particles give relatively larger con- 
tributions to the entropy density than they do to the 
susceptibility (see Eqs. (g) and @). 

In contrast, in a non- interacting QGP with 2 
massless flavors and N c colors, Rqgp—\N c ^/ (37 x 
47r 2 /90)=0.034. Perturbative calculations to higher or- 
der in the strong coupling constant can be found, e.g., 
in Ref. |18|. The use of nonperturbative results from 
lattice gauge simulations Jl9| to extract information on 
the electromagnetic susceptibility has been alluded to, 
e.g., in Ref. Q. As expected, the fluctuations in a QGP 
at temperatures T 3> T c seem to be rather close to the 
naive perturbative result. Close to T c , the situation is 
less clear. 



We turn now to estimates of charge fluctuations as- 
suming thermal equilibrium (this assumption is relaxed 
later in Sec. IV). A meaningful quantity accessible to ex- 



periments is the magnitude of the hadronic fluctuations 
normalized to the entropy S in a suitable subvolume Vs 
of the system (realized, e.g., via a restricted rapidity in- 
terval of data): 



R 



(AQ) 2 (AQ) 2 /y 3 



S 



s/v 3 



(7) 



where s is the entropy-density. In the ideal gas approxi- 



mation, 



K>) a = (0 /(<*) N/Nl 



(8) 



and 



d 3 k 
(2tt) 3 



± f(uji)lnf(uji) 

+ (l T /N))m(lT/K>] ) (9) 



where the summations run over all particle species i; the 
upper (lower) sign refers to fermions (bosons), df 1 = 
^s P m J2j 3 (ql 3 ) 2 denotes the (squared) charge degener- 
acy of a given isospin-multiplet, di is the corresponding 
total degeneracy, and uif = mf + k 2 . The thermal dis- 
tribution functions / = f B or f F correspond to bosons 
or fermions, respectively. It is important to note that, 
despite their simplicity, Eqs. (§) and (||) take account of 
interactions between the hadrons to the extent that the 
thermodynamics of an interacting system of hadrons can 
be approximated by that of a mixture of ideal gases of 
"elementary" particles and their resonances (cf. U and 
references therein). 

For a massless pion gas, 7?° = 0.253. With the 
empirical vacuum pion mass, Rjf (T c hem) — 0.18 at a 
typical chemical freeze-out temperature T c hem = 180 



D. Average Phase Space Densities 

Liouville's theorem guarantees that the final state av- 
erage phase space density of a system gets frozen at a 
certain time and gives a measure of the dynamics in the 
prior interacting stage. As observed in Ref. SJ, interfer- 
ometric (Hanbury-Brown and Twiss) studies of particle 
correlations provide a measure of the average phase space 
density of the particles in the final state. We add here 
that such knowledge also provides for a determination of 
the Debye mass. This is easily seen by recasting Eq. (^) 
for the squared Debye mass to read as 



o \2 



(m D ) 



rich 
T 



[1 T (f)ch] , or 



(f)ch = T 



T{m D f 
rich 



- 1 



(10) 
(11) 



where n c h is the total number density of charged parti- 
cles and (f)ch is the average phase space density of all 
charged particles. To the extent that n c h,T, and (f) c h 
can be determined from experiments, at least at the di- 
lute freeze-out stage of a heavy-ion collision, an estimate 
of the Debye mass is rendered possible from Eq. (|lC|). 
Equation (|ll|) is also useful, insofar as it provides a con- 
sistency check on models that predict the Debye mass 
from first principles p(| . 

For a single species of bosons or fermions, the average 
phase space density may be expressed as 



d 3 k 



d 3 k 



f 



(2tt) 3 j I J (2tt) 3 

E (±) 3+1 (lir) exp[(j + l)/3/i] K 2 [(j + l)0r 
j=i v ' 

£(±Wy exp[j/3 M ] K 2 \j0m] 

7=1 ^ ' 



(12) 



where the +(— ) sign refers to bosons (fermions) and 
Ki(x) is a modified Bessel function. For temperatures 
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such that T/m <C 1, use of the asymptotic relation 
K v {x) ~ (n/2x) 1 ^ 2 exp(-ar) reduces Eq. @ to the re- 
sult 



B,F\ 



exp[/3(/i — m)] 
(271) 



(13) 



which is valid in the classical regime. 

For orientation, we show in Fig. [l] results for the av- 
erage phase space density (f B ) c h of charged pions and 
kaons for selected values of their respective chemical po- 
tentials as functions of temperature. Although results 
for temperatures up to 200 MeV are shown in this figure 
for illustrative purposes, we should focus our attention 
on temperatures below T c , the temperature at which the 
transition to a QGP phase occurs. For reference, the dot- 
ted line in the upper panel shows the result (f B ) c h — 0.37 
for massless pions at zero chemical potential (in contrast, 
(f F )ch — 0.088 for massless fermions at zero chemical 
potential). All other curves show results with vacuum 
values of masses. The dashed curves in both panels are 
for zero chemical potentials. Evident features from the 
results of this figure are: (1) Increasing masses of bosons 
reduce (f) c hl (2) Positive (negative) chemical potentials 
of bosons enhance (suppress) (f B ) c h relative to its value 
at /i = 0. For chemical potentials approaching the mass 
of the particle (incipient Bosc condensation), the aver- 
age phase space occupancy begins to become increasingly 
large; and (3) The higher (lower) the temperature, the 
higher (lower) is (f B ) c h- 
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FIG. 1. Average phase space densities (f B )ch of charged 
pions (upper panel) and kaons (lower panel) as functions of 
temperature at the indicated values of their respective chem- 
ical potentials. The dotted curve in the upper panel shows 
(f B )ch — 0.37 for massless pions at zero chemical potential. 

It is straightforward to compute (f) c h in a mixture of 



hadrons through the use of Eq. ( |l2| ) by summing over the 
various species with their appropriate charge degenera- 
cies df 1 . We refer the reader to Refs. (^Jl^l, where the 
extent to which the conditions at the freeze-out stage of 
the collision can be delimited are discussed. 



E. Electromagnetic Correlations 

An intimate connection exists between the fluctuations 
of the charged constituents in the system and the emis- 
sion of electromagnetic radiation. The e.m. polarization 
tensor figuring in Eq. (|^) directly governs the production 
rates of dileptons and photons. This is important, since 
electromagnetic emission occurs throughout the evolu- 
tion of the system created in heavy-ion collisions and 
ranks highly among the principal ways in which the dense 
stages of evolution can be directly probed (cf. |2l[] and 
references therein). 

Beginning with dilepton emission, we note that 

the rate of production from a thermally equilibrated sys- 
tem is given by 



dR 



l+l- 



d 4 q 



-^2/ S (9o;T)i(lmn i + 2Imn r ) , (14) 



where is the usual lepton tensor, and the electromag- 
netic hadron tensor 



W™{q) 



d 4 xe iq - x 



{j; m {x) Jr(0)) c (15) 



is given directly in terms of the imaginary part of the 
retarded photon self-energy 



WT = -2 f B (qo;T) Im U, 



(16) 



Photon emission rates are obtained by taking the limit 
M 2 = q 2 . — q 2 — > in the above expressions. Note that, 
unlike the space-like limit involved in Eq. (^|), dilepton 
(photon) production probes the time- like regime (limit), 
as well as both the longitudinal and transverse parts of 
the electromagnetic polarization tensor. 

We turn now to elaborate on this interrelation, begin- 
ning with the simple and transparent case of an ideal gas 
of massless constituents (either pions or quarks). In this 
case, the leading temperature contribution to the one- 
loop self-energy is (see, e.g., Ref. fll8|] ) 



IL™(q ,q) = -mj 



1-* 




9o+9 


2q 




9o-9 



(-M 2 ) 



(17) 



where g?q = 1 {N c J2 q e q) f° r pions (quarks). Note that 
n°°(9o>9) acquires an imaginary part only for space- like 
M 2 < 0, due to Landau damping. The dilepton produc- 
tion rate, which stems from the time-like region charac- 
terizing decay contributions (emissivities) , takes the form 
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dRi 



+i- 



a 2 d Q 



d 4 q 



9= o a 2 d,Q 



T S M 2 



rim, 



M 2 
48^ 



1± 



2T, 
— In 

q 

2/(|) 



(18) 



where d,Q = 1 for pions or N 2 N C J2 q e q f° r quarks (iV s = 
2), and x± = exp[-(g ±q)/2T]. 

For illustrative purposes, the q — > limit has been 
taken in the rightmost expression. Note that the overall 
Bose factor is not part of the photon self-energy so that 
the leading temperature contribution to dilepton produc- 
tion is due to the imaginary part of the vacuum piece. In 
the q — ► limit, the explicit form of the 00-componcnt 
for time-like M 2 > is 



ImU 00 (q ,q) = -^- 

1Z7T 



1± 



v(f) 



(19) 



A formal connection between the space-like and time-like 
regimes can, in principle, be constructed via a dispersion 
relation. For the case of interest here, the analytic prop- 
erties of the retarded photon self-energy imply that 



RelP o (q = 0,q) = -J ^ 



duj 2 Imn 00 (w,g) 



-{ui 2 



(L 2 Imn 00 (w,g) f duj 2 lmU m (uj,q) 



7T -(CU 2 ) 



, (20) 



where the second equality explicitly shows the decompo- 
sition into space-like and time-like contributions. Using 
the imaginary parts from Eqs. ( p"7j ) and (jll]), we find 
that Ren 00 (qo = 0, q), and hence the Debye mass, is sat- 
urated by the space-like part of the dispersion integral, 
the second term in Eq. (|20|) giving a vanishing contribu- 
tion. Thus, for a thermally equilibrated system at rest, 
a model-independent relation between fluctuations and 
dilepton production cannot be achieved. However, within 
a given approach for calculating electromagnetic emissiv- 
ities, the Debye mass and hence charge fluctuations are 
simultaneously determined. 



F. Effects of Interactions on Electromagnetic 
Correlations 

The interrelations between fluctuations, susceptibili- 
ties, and the photon polarization function described in 
the previous sections set the stage to address the effects of 
interactions between the constituents in a thermal equi- 
librium system. In a hadronic gas at a finite temperature 
and baryon chemical potential, various approaches to cal- 
culate W^" have been pursued 21 26|. In the following, 
we employ two examples related to dilepton and photon 
emissivities to analyze the impact of interactions on fluc- 
tuations. 



1. Vector Dominance Model 

A number of hadronic models used to study medium 
modifications invoke the phenomenologically success- 
ful Vector Dominance Model (VDM) at finite tem- 
peratures and densities using effective hadronic La- 
grangians ^||,|3). In VDM, the photon couples to 
charged hadrons exclusively through the vector mesons 
p, uj and <p. In the isovector p-channel, which constitutes 
the dominant part of the photon self-energy, the latter 
then takes the form 



U^(q ,q;fx B ,T) 



(0)n4 



9l 



■D^(q Q ,q;p Bl T). (21) 



With the help of standard projection operators P£ v and 
Py" [ fl8f , the in-medium p propagator can be decomposed 
into its longitudinal (L) and transverse (T) parts: 



D^(q Q ,q;p B ,T) = D L p Pr 



D T p P^ 



D 



L,T _ 



M 2 



(m^) 2 - E L P ' T 



(22) 



Above, g p denotes the universal VDM coupling constant, 

m{,°' is the bare p mass, and Sp :T stands for all possible 
proper self-energy insertions. 

Let us first focus on the p — > irir loop contribution to 
S p , which involves the strong coupling of the p. The re- 
summation of this loop to all orders via the propagator 
accounts for interactions, albeit from a select class of di- 
agrams. To extract the electromagnetic Debye mass, we 
calculate the prnr loop at finite temperature |2^] , take the 
space-like limit and, after subtracting the vacuum piece, 
obtain 



m n = -- 



7 ( °V 

9 ' km [D™(0,q)-D 00 

q^,0 



9 



o ^2 



p,vac 



(0,0)] 



(4 0) ) 2 



(mp 0) ) 2 + m 2 D 



(23) 



with (m° D ) 2 given by Eq. (||). The quantity rh 2 D ~ 
9p ( m r>) 2 represents the "strong" (isovector) Debye mass 
induced by the pTnr loop. The rescattering of pions 
through the formation of the p resonance causes a re- 
duction of the Debye mass; in turn, the hadronic charge 
fluctuations are reduced. This result is shown by the 
solid line in Fig. where at temperatures close to T c h em 
a reduction of ~ 25% is apparent. If one accounts for 
effective pion-number conservation (which is implied by 
chemical freezeout in a heavy- ion collision), the resum- 
mation effect is more pronounced at lower temperatures. 
At SPS energies, pion-number conservation entails the 
build-up of a finite pion chemical potential which can 
reach values p^ ~ 70 MeV for thermal freezeout at 
T fo ~ 120 MeV (25). This effect is illustrated by the 
dashed curve in Fig. H, which shows that even around 
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thermal freezeout a reduction of ~ 20% of the fluctua- 
tion content persists. 

Contributions from heavier states in a hot mesonic 
environment are also present, but are less significant. 
The KK loop enters predominantly into the (f> meson 
component of VDM (OZI rule), and the pp loop is 
kinematically suppressed. For typical T c hem, we find 
(mD, w /mD,ii) 2 -l/'l, which when resummed results in 
additional reduction, but only at at the few percent level. 
Potentially more important are "off-diagonal" interaction 
terms, most notably thermal 7ra 1 (1260) loops which are 
induced by resonant p scattering off thermal pions, or, 
if abundant, baryonic excitations. Various interaction 
vertices have been employed in the literature, see, e.g., 
Rcf. (27]] for a recent survey. It turns out that the fre- 
quently employed p-meson tensor (p^) couplings, which 
have the advantage of individually maintaining gauge in- 
variance, lead to a vanishing contribution to the electro- 
magnetic Debye mass. However, other choices, which do 
not vanish in the space-like screening limit are possible. 
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FIG. 2. VDM resummation correction to AQ 2 represented 
by the ratio of the squared electromagnetic Debye mass that 
includes effects of resummed pirn loops (according to Eq. (p3[)) 
to that from the lowest-order result in strong interactions. 



2. Chiral Reduction Scheme 

The framework of Ref. |24| allows chiral master formu- 
lae to be used to relate to experimentally extractable 
spectral functions in a model independent way. Using the 
analytical properties of Woo i n the complex u-plane, it 
follows that, in thermal equilibrium, 



Woo (<".?) 



2lmIl o6 ab ( UJ 'i) 



eP\w\ + 1 

1 00 



a/3" 



1 



(24) 



where the superscript F (R) explicitly indicates the Feyn- 
man (retarded) charge-charge correlator. In this scheme, 



the vacuum part H 0F is common to a hadronic or par- 
tonic gas and is related to experimentally measured spec- 
tral functions. One has 



n 



00 



(M ) = i / d*x<?** (0\T* VS(x) V o b (0)|0) c , (25) 



which is given by the longitudinal part of the isovector 
correlation function, 



tS ab (g 00 M 2 ~q 2 )n^(M 2 ) 



(26) 



In a hadronic gas, the corrections to the leading spectral 
density follow from a virial expansion. Specifically, 



nlF.ab 1 \ 
00 ( w ><?) 



E 

A 



d 3 k r(uj k ) 



(2tt)3 2u k 
xi I d A xe iq ' x (n(k)\T* V£(x) V*(0)\n(k)) t 



+..., (27) 

where lo\ = m 2 . + k 2 . The ellipsis refer to corrections of 
higher order in pion- and nucleon-densities. The reduc- 
tion of the forward Compton amplitude in Eq. (^) in- 
volves both the vector and axial-vector correlation func- 
tions. One finds 



ImIl^' ab (co,q)^S a 



d 3 k r(tu k ) 



(2tt)3 2/2 UJ k 
-Aq 2 lmll v (q) 

+2{k + q) 2 Im n A (k + q) + (q -> - 
+4(2wfe + u) {uj k q 2 - u k ■ q) 

xlmIl v (q)Bje(A R (k + q)+q^> 



■q) 



(28) 



where / w is the pion decay constant and A#(fc) is the re- 
tarded pion propagator. In obtaining this result, only the 
dominant contributions to the Compton amplitude were 
retained. We note that the current algebra result follows 
from setting k = 0, as a result of which the quantity in 
brackets in Eq. (^8|) takes the form 

(29) 



Aq 2 (ImU A (q) - ImU v (q)) 



This relation illustrates the role of chiral symmetry 
restoration at next-to-leading order in the hadronic (pio- 
nic) density. Higher order corrections, not discussed here, 
are calculable in a similar way and in close analogy with 
calculations of photon and dilcpton cmissivitics. 



III. CONNECTIONS TO HEAVY-ION 
EXPERIMENTS 

Since fluctuations at the dilute freeze-out stage of a 
heavy-ion collision have been amply discussed in the lit- 
erature, we focus here on how one might gain knowledge 



G 



about fluctuations during the dense stages of the colli- 
sion. Direct photons and dileptons, due to their negligi- 
ble final state interaction in a heavy-ion collision, have 
long been recognized as sensitive probes of the physi- 
cal state of strongly interacting matter. We are thus 
led to the following natural question: To the extent that 
electromagnetic observables at the CERN-SPS have been 
adequately described through a combination of hadronic 
and partonic approaches |2|j2^,^] (see |2l| for a review), 
what are the implications of dilepton and photon spectra 
to be measured at RHIC and LHC on charge fluctuations 
of hadrons and partons during the dense stages of evolu- 
tion? The answer to this question hinges critically on the 
outcome of upcoming measurements and on the ability 
of theory to model the complex evolution of the system 
created in a heavy-ion collision. In the following section, 
we provide a suitable theoretical framework that can be 
systematically improved. 



A. Longitudinally Expanding System in Local 
Thermal Equilibrium 

In a heavy-ion collision, isospin fluctuations occur 
in an expanding environment that initially consists of 
partons which hadronize towards the end of evolution. 
In the Bjorken expansion scenario with boost invari- 
ance [^9| , the isospin charge per unit pseudo-rapidity can 
be schematically defined as dQ a /dr]. Boost invariance 
and charge conservation imply that 



"7+ 



d(dQ a / drf) 
dr 



0. 



(30) 



This conservation law assures that the mean fraction of 
isospin is the same in all frames, and constant through- 
out the proper time evolution. If this were to hold for 
the variance, then we would naively conclude that the 
rest frame results developed in the preceding section will 
trivially extend to the rapidity-binned fluctuations in a 
heavy-ion collision. Closer inspection reveals that this is 
not the case, as will be shown below. 

Consider the thermal fluctuations in an expanding 
fluid-element that is characterized by its pseudo-rapidity 
rj = tanh~ (t/z) and proper time r = \Jt 2 — z 2 . Follow- 
ing the standard generalization of a conserved current to 
a moving frame, we define the isospin charge operator 



V" 



(31) 



in a fixed pseudo-rapidity interval of extent A centered 
around 77. Above, represents a 3-dimensional hyper- 
surface in space-time which is orthogonal to the local 
fluid 4-velocity = 7(1, (3). From here on, we restrict 
ourselves to purely longitudinal expansion for the sake 
of simplicity. With da^ = dxj_(dz,Q±,dt), the charge 
becomes 



Q1(^V)= / dx ± / rdr] uZV£{Tch.T),x ± ,TshT)) , (32) 



where rj± = r\ ± A/2. The quantity tifj 1 = (ch??, 0_l, shrj) 
is equivalent to the longitudinal Lorentz boost operator 
Ag^, as it should. The binned isospin fluctuations can 
then be expressed as 



v+ v+ 



(CaM) QaM)>c = J dx x dx' ± J J rd^rdr/ 



x u(f u'l; {Vlix^xMix'^x'^c, (33) 

where xii = r(chry, shr;). The averaging is carried in a 
hadronic or a partonic gas in equilibrium at temperature 
T(t) = 1//3(t) and pertinent chemical potentials. (We 
will return to discuss some specific issues related to non- 
equilibrium situations in Sec. IV.) Using translational 



invariance in the transverse direction along with a Fourier 
transform of the unordered charge correlator, 



(34) 



W$(u,q) = J d 4 x e** (V °(x) V o b (0)) , 

Equation ( |33| ) can be cast in the form 

du) dq z 



(Ql(r,v) Qa(t,v))c = V_ 
x Wf(u>/r,0 ± ,q z /r) 



(2n) 2 



J 2 -q 2 z 



(35) 



Here, we have introduced dimensionless variables lj = tlo 
and q — Tq and exploited the relations 

W 0Q = ^W L , W Qz = ^W L , & ndW zz = ^W L , 

as required by current conservation. The kinemati- 
cal "formfactors" in Eq. ([55|) arise from restricting the 
isospin charge to a bin of width A in pseudo-rapidity 
and are given by 



Fa Ft 



2 sin [ sh ( — ) [to shrj — q z chrj] 



(36) 



The occurrence of solely the longitudinal part Wl is a 
consequence of the underlying assumption of purely lon- 
gitudinal fluid flow and current conservation. The pres- 
ence of transverse flow would naturally involve the trans- 
verse part Wt as well. 

We recall that for dilepton emission, the combination 
2 Wt + Wl contributes, whereas in photon emission only 
Wt is needed. At large proper time r, the correlation 
function is probed around to « q w 0. In general, the fluc- 
tuations in pseudorapidity bins represented by Eq. ( |35| ) 
are sensitive to both space-like and time-like physics. As 
noted in the previous section, the space-like contribu- 
tion is driven by screening and Landau damping, while 
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time-like physics governs decays or production. The lat- 
ter are at the origin of photon and dilepton emissivities. 
To lowest order in temperature, the imaginary part of 
the relevant polarization is empirically accessible through 
e + e~ annihilation for M 2 > 0. Since thermal fluctua- 
tions preserve isospin, the integrand in Eq. ( |35| ) (modulo 
F) is a consequence of the fluctuation dissipation theo- 
rem (cf. Eq. (|24)). 

Employing boost invariance, we can focus on a window 
of width A centered around mid-rapidity rj — 0. Insert- 
ing Eq. ( |2^ ) into Eq. (|33|), we obtain (omitting isospin 
indices) 



= (-2) [2sh(A/2)r 
L,q/r) 



1 



duj dq 

sin [gsh(A/2)]V 
gsh(A/2) J 



(37) 



This relation is one of the principal results of this work. It 
explicitly demonstrates that fluctuations are sensitive not 
only to screening, but also to emissivities. Such a feature 
is in contrast to the case of the rest-frame fluctuations 
discussed in Sec. [FJ, and stems from the general properties 
of boost transformations of charge-charge correlators. 

Equation (^) also shows that fluctuations are pro- 
portional to the "size" of the pseudorapidity interval 
through sh 2 (A/2), vanishing for A = as they should. 
The dominant contribution arises from long wavelengths 
A ~ 2n/q > 27rrsh(A/2) for which sin 2 x/x 2 fts 1, 
and small frequencies / = lo/2ttt ~ l/2nTp, where 
Tp is a typical lifetime of the expanding system. The 
limitation on frequency is enforced by the requirement 
that for large times, one has to recover the equilibrium 
screening limit that was discussed in the previous sec- 
tion. That this is indeed the case can be seen as fol- 
lows: For large T,the oscillating factor on the right- 



hand-side of Eq. d37|) collapses the q integration to its 
saddle point at q ~ 0. Furthermore, the denominator 
becomes e l3u) / T — 1 ~ /3lj/t. Using this, one recovers the 
dispersion relation, Eq. (|20|), for the Debye mass. Note 
also that the factor z — 2rsh(A/2) combines with V± to 
yield the spatial 3- volume V3. 

Figure displays results of the quantity R = 
(Qa( t ))/^S s ( t ) f rom numerical evaluations of Eq. d37j ) 
and Eq. (^|) for a model system comprised of an ideal 
QGP (Nf = 2) in the partonic phase and massless pions 
in the hadronic phase. The time evolution of temperature 
has been taken to result from a boost-invariant longitu- 
dinal expansion |29fl in which 



To (ro/r) 1 ^ 
T(t) = { T c 



r <r<r 1 (QGP) 

n < t < t 2 (Mixed) (38) 

T 2 < T (HG) . 



The numerical values of the time constants, tq = 1 fm/c, 
Ti = 2 fm/c, r 2 = 5 fm/c and the critical temperature, 
T c = 170 MeV, were chosen to generate conditions that 
resemble those at SPS energies (e.g., T — T^ti/to) 1 / 3 ~ 



214 MeV). In the mixed phase, the fluctuation content 
R = JRqgp + (1 — J)Rtt is a weighted sum with the 
appropriate volume fractions / and (1 — /) of the QGP 
and hadronic phases, respectively. 



u 

13 

o 

'l/j 

c 

a 




equilib. limit 

— space-like (free) 

— time-like (free) 

— 7T7T with F vac 

TC 

... „ med 
7171 With t 
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FIG. 3. Proper-time evolution of charge fluctuations per 
unit 3- volume, Eq. _d37| ) , normalized to the appropriate en- 
tropy density, Eq. (BhTfor an ideal QGP (r < 2 fm/c) and a 
massless pion gas (r > 5 fm/c) in a locally thermalized, lon- 
gitudinally expanding fluid. The upper (lower) panel shows 
results for a pseudorapidity window of size A = 1 (0.5). The 
solid (long-dashed) curves correspond to contributions from 
the space-like (time-like) part of the charge-charge correla- 
tor. The dash-dotted and double-dash-dotted lines are re- 
sults for the time-like hadronic contribution wherein nn cor- 
relations have been included through the free pion-e.m. form- 
factor F,T C (M 2 ) and an in-medium F™ ed (M 2 ) (with twice 
larger p- width), respectively. The dotted lines indicate the 
static equilibrium limit, Tm 2 D /s. 

In the non-interacting case (employing for Im Wp° 
Eqs. (|l7|) and ( |l9| ) generalized to finite q), fluctuations 
are clearly dominated by the space-like contribution of 
the correlator (solid curves in Fig. |3j), especially at small 
proper-times in the partonic phase. The latter fact is 
due to the fermionic character of the cha rge c arriers (see 



Eq. (19)). However, as shown in Sec. |II F| , the inclu- 
sion of hadronic interactions tends to suppress the space- 
like contribution (not shown in the figure) , but enhances 
the time-like one, as is evident from the dash-dotted and 
doublc-dash-dotted curves in the upper panel. This be- 
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havior becomes more pronounced if the size of the rapid- 
ity bin is reduced (see the lower panel of Fig. |j). 

In general, A should be large enough to enhance the 
volume-to-surface effect, but small enough to avoid the 
constraint of overall charge conservation in the nuclear 
system. 

We emphasize that the evolution of fluctuations shown 
in Fig. H refers to a highly idealized model, especially in 



the hadronic phase. As stressed in Sec. IIC, the en- 



tropy density of a more realistic hadronic resonance gas 
will be several times larger than that of a massless pion 
gas ||, entailing a reduction in R. To illustrate its ef- 
fect, we show in Fig. [| the results of a calculation in 
which the pure hadronic phase begins with an entropy 
density s ~ 5 fm~ 3 , a value representative of a hadronic 
resonance gas at T c = 170 MeV. Finite hadron masses 
also imply that sjjg falls off faster than T 3 with decreas- 
ing temperature, which has been simulated by employing 
s hg( t ) = s hg(T c ) ^/t) 1 - 35 , indicated by thermal fire- 
ball calculations 123]. Consequently, the static limit of 
the ratio R increases with falling temperature in the ex- 
panding resonance gas for r > 5 fm/c (see the dotted 
lines in Fig. ||). To minimally account for the effect of 
additional states in the numerator of R (this contains the 
magnitude of fluctuations) , we have also included contri- 
butions from kaons to the hadronic e.m. correlator (in 
the massless limit). This amounts to a simple increase of 
(Qa( t ))/^3 by a factor of 2 compared with the result for 
massless pions. The presence of massless strange quarks 
in the QGP hardly affects the equilibrium value of Rqgp, 
reducing it from 0.034 (for N f = 2) to 0.032 (for N f = 3). 
The relative role of space- and time-like contributions to 
R does not change as compared to Fig. ||. 

More realistic hadronic correlators, as were found nec- 
essary to explain the dilepton data from heavy-ion col- 
lisions at SPS energies, are characterized by a substan- 
tial reshaping of the pion e.m. form factor, in particular 
through an appreciable accumulation of strength at low 
invariant masses. We anticipate that the use of such cor- 
relators in Eq. (|37]) will enhance the role of the time-like 
region. Results for the simplified cases considered here 
are thus to be regarded as suggestive. They clearly indi- 
cate, however, that calculations of a more realistic evolu- 
tion of R employing improvements in (i) the equation of 
state, (ii) space-time evolution, and (iii) appropriate elec- 
tromagnetic correlators are worthwhile. Developments 
in this regard are in progress and will be reported sepa- 
rately. 




4 6 8 10 

T [fm/c] 

FIG. 4. Same as Fig. |^ but with the following improve- 
ments: (i) In the hadronic phase, the entropy density has been 
calculated for a hadronic resonance gas, (ii) The e.m. correla- 
tor includes contributions from K + K~ states in addition to 
tt + tv~ states (both are calculated in the massless limit); and 
(iii) In the QGP, massless strange quarks have been included 
in both the entropy density and e.m. correlator. 



B. Relation to Other Works 

Our main result from the previous section, Eq. ( |37|) , 
bears some resemblance to the results obtained in Rcf 7M , 
in particular to Eq. (18) therein, which, in the notation 
of Ref . H , reads 



«AQ)*> = /| 



dk sin (kAy/2) 



{k/2f 



Xk 



:[l-cxp{-fc 2 (Ay diff ) 2 }] 



(39) 



This relation describes the relaxation of the moments 
of charge fluctuations in one-dimensional rapidity space 
through the purely space-like density-density correlation 
function Xk and the mean diffusion, A^diff, of a particle 
during time r. The (sinx/x) 2 factor above is equiva- 



lent to a similar factor in Eq. (37) and arises as a conse- 
quence of boost invariant longitudinal expansion without 
transverse flow. The time-dependent exponential factor 
in Eq. ( |39| ) results from non-equilibrium effects, which 
were treated using a Gaussian solution of a Fokkcr-Planck 







equation S. Non-equilibrium effects are not accounted 
for in Eq. ( p7|) (see, however, our discussion of such effects 
in the following section). Our derivation of Eq. (|3^), on 
the other hand, includes the full 4-dimcnsionally trans- 
verse structure of the conserved isospin current operator, 
evaluated within a locally equilibrated relativistic fluid. 
This induces additional quantum effects associated with 
the finite time scale of the expansion thereby involving 
the time-like region of the correlation function. These 
features are not captured by the treatment of Ref . || , in 
arriving at Eq. ([59). 

It is worthwhile to point out that, during the dense 
stages of evolution in a heavy-ion collision, the assump- 
tion of local thermal equilibrium among the constituents 
of matter has been gainfully employed to understand 
photon and dilepton emission. Calculations in the CERN 
SPS regime have shown that a fair agreement exists be- 
tween theory and experiment |25|]26] , |2"8| ] . In the context 
of charge fluctuations the equilibrium assumption is sup- 
ported by the estimates of Ayaiff presented in Ref. 
and also by the more quantitative analysis of Ref. [gj, 
where the diffusion equation ( |39| ) was solved with inputs 
from microscopic transport calculations || . It was shown 
there that, due to frequent rescattering in the hadronic 
phase of a heavy-ion collision, the diffusion of charge in 
rapidity space is effective in keeping fluctuations close to 
its equilibrium value (within ~ 20%) until freeze-out at 
both SPS and RHIC ||. This is consistent with prelim- 
inary experimental data, which point to a dilute hadron 
gas in the vicinity of thermal freeze-out. 

Notwithstanding these observations, we study in the 
following section a somewhat different source of non- 
equilibrium effects, namely, fluctuations in the temper- 
ature. This was suggested some years ago as a means 
of extracting the heat capacity of the hadronic gas close 
to freeze-out Jll]]. Here we explore connections between 
temperature and charge fluctuations, and examine the 
consequences of the presence of non-equilibrium effects 
in heavy-ion collisions. 



IV. DEVIATIONS FROM LOCAL THERMAL 
EQUILIBRIUM 

Non-thermal fluctuations are generally caused by com- 
plex microscopic and dissipative phenomena. A micro- 
scopic or QCD-based description of such phenomena is 
beyond the scope of this paper. Instead, we use a macro- 
scopic approach that is based on a stochastic or maxi- 
mum noise assumption. We adduce arguments to show 
that a non-thermal or early-stage isospin fluctuation can 
survive the expansion only if it is accompanied by a large 
temperature fluctuation at freeze-out. Since temperature 
fluctuations reflect the heat capacity of the underlying 
phase, we consider two different possibilities by treat- 
ing non-equilibrium effects stochastically through addi- 
tive (Gaussian) or multiplicative (power-law) heat flow. 



A. Gaussian Fluctuations 

Let T(t) be the temperature at a given proper time r, 
and T* = 1//?* the equilibrium (mean) temperature. Let 
T(t) relax stochastically due to heat flow. In a phase 
in which the heat capacity C depends strongly on T*, 
thermodynamics implies that the variance of the tem- 
perature distribution is bounded. The stochastic evolu- 
tion is "Brownian" , and the relaxation of T(r) is gov- 
erned through additive white noise. The corresponding 
Langevin equation is of the form 

^ + - (T(r)-e(r)-T,)=0, (40) 
or r* 

where r + is the temperature relaxation time, and £(r) is 
a white noise, i.e., uncorrelated in time and with zero 
mean. In this case, 

«£(tK(t')» = 2D*(t-t'). (41) 

The average in Eq. ( pll| ) is performed over stochastic noise 
distributions £, the strength of which is characterized 
by the variance D. The distribution of temperatures 
PJt, T) associated with the additive Langevin equation 
([i"o| ) obeys a Fokker-Planck equation . Specifically, 

^T) = -—\—— P(r,D 



1 8 2 (2D 



(42) 



with a drift rate vt = (T*—T) / t* and a constant diffusion 
rate u^iS — 2D/t 2 . The general solution is an evolving 
Gaussian 



P(t,T) = 



1 



-(T-(T(T))) 2 /2* 2 (r) 



2ir cr 2 (r) 

with a stochastic mean 

(T(t)-T,) = {T(0)-T,)c- t / t *. 



(43) 



(44) 

Above, single angular brackets denote an average with 
respect to the temperature distribution. The variance is 
given by 



<r 2 M 



£ + „*(0) 

T* \ T, 



£.\ „-2t/t. 



(45) 



which reduces to <t 2 (t — > oo) = D/t* in the stationary 
case, reflecting the thermodynamic limit. It can thus be 
related to the heat capacity, C = T 2 /a 2 , which depends 
strongly on T„ . The stationary temperature distribution 
takes the form 



P s (T) = P(oc,T) 



C 



2ttT 2 



-C(T-T,) 2 /2T? 



(46) 



The standard Gaussian nature of the fluctuations attests 
to the additive character of white noise. Obviously, a 
large heat capacity strongly limits the character of fluctu- 
ations through exponentially small tails in the stationary 
distribution, Eq. @. 
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B. Power-Law Fluctuations 



V. SUMMARY AND OUTLOOK 



Here we explore the consequences of a more speculative 
scenario, that is, a phase in which the heat capacity C 
is nearly constant so that the equilibrium variance grows 
with the mean temperature squared, a 1 = T^/C. Such a 
behavior can be constructed by starting from a Langevin 
equation in which the noise is introduced multiplicatively 
using the Ito-prescription |[l) . In our case, this is realized 
by setting 

^ + ~\ [(l + e(r))T(r)-T,]=0, (47) 

where, as before, £(t) is a white noise with zero mean 
and variance 2D. Arguments similar to the ones used 
in the previous section yield a Fokker-Planck equation 
analogous to Eq. p^ ) with the same drift rate, vt — 
(T* — T)/t», but a non-constant diffusion rate, v^is = 
2DT' 1 1 't„ 2 , which is the key difference from the result in 
Eq. ((iCf). In the steady state, the distribution of temper- 
atures is now given by 



Pscn = 



l+C 



(1 + C) 1+I 

r(i + c) \t 



T \ 3+C e -(l + C)T«/T 



I 



(48) 



where the corresponding variance, a 2 — T 2 /(t*/D — 1), 
allows us to replace t,/D by 1 + C. Note that this is 
the desired behavior for a constant heat capacity. Larger 
temperatures are only power-law suppressed as a direct 
consequence of the multiplicative character of the noise. 
If we further assume that the temperature relaxation 
time Tit is small in comparison to the underlying flow 
time in the evolution, we may average the temperature 
fluctuations around the instantaneous mean = 1//3* 
using Eq. pH). For the leading temperature contribu- 
tion to the electromagnetic correlator, e.g., Eq. (J27n, this 
amounts to 



l 



■j/T _ l 



n=l \ 1 + n T, (l+C)/ 



(2+C) 



(49) 



Again, we see that temperature fluctuations turn the ex- 
ponential suppression of high energies (implicit in equi- 
librium Bose and Fermi distribution functions) into a 
power-law suppression. The larger the heat capacity C 
of the phase considered, be it hadronic or some sort of 
mixed phase, the stronger will be the power-law suppres- 
sion, and hence the closer the temperature distribution 
to a Gaussian, as expected. These results can be readily 
adapted to a dynamical description, e.g., a hydrodynam- 
ical simulation, which will be reported in future work. 



Our purpose in this paper has been to highlight the in- 
terdependence of fluctuations in isospin charge and elec- 
tromagnetic emissivities in the strongly interacting mat- 
ter encountered during the dense stages of a heavy-ion 
collision. We have established such connections both in 
an infinite system at rest, as well as in a relativistically 
expanding fluid to mimic some of the dynamics present 
in high-energy nuclear collisions. 

In the static case, we have emphasized that well-known 
relations between the fluctuations of the charged con- 
stituents of matter and its electric screening mass im- 
ply an intimate connection to the emission of dilep- 
tons/photons. Both of these can be calculated from the 
retarded current-current correlation function, or the pho- 
ton self-energy, in hot and dense partonic or hadronic 
matter. Consequently, analyses of dilepton and photon 
observables that are based on the e.m. correlation func- 
tion can be systematically applied to study medium ef- 
fects on fluctuations during the dense stages of a heavy- 
ion collision. Two hadronic approaches were used to ex- 
emplify that close to the expected phase boundary, hot 
and dense hadronic matter might exhibit fluctuations of 
similar magnitude as an equilibrated partonic medium 
(QGP). Clearly, further work is required to improve these 
estimates. 

In the case of an expanding fluid, we derived an expres- 
sion for the fluctuations in a locally equilibrated fluid ele- 
ment as a function of proper time and over a finite range 
in pseudo-rapidity. Isospin conservation, which governs 
the 4-dimensionally transverse Lorentz-structure of the 
correlation function, implies that fluctuations depend on 
both transverse and longitudinal components. Further- 
more, the finite time scale ~ r/ G dictated by the strong 
interaction dynamics in a nuclear collision implies that 
both space-like and time-like processes are involved. At 
a given freeze-out proper time r/ and within a rapid- 
ity window of width A, the thermal isospin fluctuations 
are dominated by constituents with typical frequencies 
/ ps 27t/t/ and wavelengths A » 7tt/ sh(A/2). 

Our analysis of non-equilibrium effects associated with 
possible temperature fluctuations indicates that a large 
isospin fluctuation at the early partonic stage can survive 
the course of a relativistic heavy-ion collision only if a 
large temperature fluctuation occurs simultaneously at 
freeze-out in the same event. Given a large and strongly 
temperature-dependent heat capacity of a hadronic gas, 
this scenario seems unlikely. 

Prospects for the future include: (1) A determination 
of the Debye mass of the system at the dilute freeze- 
out stage of a heavy-ion collision; this could be accom- 
plished solely through measured results of particle multi- 
plicities and particle correlations using Hanbury-Brown 
and Twiss analyses. An intriguing question in this con- 
nection is whether or not the Debye mass at freeze-out 
remains constant from AGS through LHC energies; and 
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(2) A delineation of the role of charge fluctuations during 
the dense stages of the collision through a study of elec- 
tromagnetic emissivities; this is a more challenging task 
inasmuch as an understanding of the measured dilepton 
and photon yields will necessarily require a detailed the- 
oretical analysis to unravel the phase structure of the 
rapidly evolving produced matter. 

Our ongoing investigations include a dynamical simu- 
lation of the collision, in which the effects of strong in- 
teractions are studied through the equation of state and 
the effects of fluctuations in dense matter are studied 
through its electromagnetic radiation. 
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